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Abstract. We present ESO/VLT spectra in the 2.9-4.1 ^m range for a large sample of infrared stars in the Large 
Magellanic Cloud (LMC), selected on the basis of MSX and 2MASS colours to be extremely dust-enshrouded 
AGB star candidates. Out of 30 targets, 28 are positively identified as carbon stars, significantly adding to the 
known population of optically invisible carbon stars in the LMC. We also present spectra for six IR-bright stars 
in or near three clusters in the LMC, identifying four of them as carbon stars and two as oxygen-rich supergiants. 
We analyse the molecular bands of C2H2 at 3.1 and 3.8 fim, HCN at 3.57 pim, and sharp absorption features in 
the 3.70-3.78 /xm region that we attribute to C2H2. There is evidence for a generally high abundance of C2H2 
in LMC carbon stars, suggestive of high carbon-to-oxygen abundance ratios at the low metallicity in the LMC. 
The low initial metallicity is also likely to have resulted in less abundant HCN and CS. The sample of IR carbon 
stars exhibits a range in C2H2:HCN abundance ratio. We do not find strong correlations between the properties 
of the molecular atmosphere and circumstellar dust envelope, but the observed differences in the strengths and 
shapes of the absorption bands can be explained by differences in excitation temperature. High mass-loss rates 
and strong pulsation would then be seen to be associated with a large scale height of the molecular atmosphere. 

Key words. Stars: carbon - Stars: AGB and post-AGB - circumstellar matter - Stars: mass-loss - Magellanic 
Clouds - Infrared: stars 



1. Introduction 

One of the main contributors of elements such as carbon 
and nitrogen, and arguably the most important "factory" 
of cosmic dust, Asymptotic Giant Branch (AGB) stars 
represent the final evolutionary stages of intermediate- 
mass stars (MzAMS ^ 1 to ~ 8 Mq), when they lose up 
to 80 per cent of their mass at rates of AI ^ 10~® to 10^^ 
M0 yr-i. Aged between t ~ 30 Myr and t > 10 Gyr, AGB 
stars chemically enrich the Universe on timescales ranging 
from less than the dynamical timescale of a galaxy to as 
long as the age of the Universe. AGB mass loss happens 
as a result of two mechanisms: (1) strong radial pulsations 
elevate the stellar atmosphere out to a distance where (2) 
the temperature is sufficiently low {T 1000 K) and the 
density is sufhciently high {n-a^^lO^^ m"^) for dust for- 
mation to occur. Radiation pressure from the luminous 
giant {L ~ 10^ Lq) drives away the dust, taking the gas 
along with it by mutual collisions. During this "super- 



wind" stage the dust-enshrouded AGB star vanishes at 
optical wavelengths but shines brightly in the infrared. 

An unsolved problem of AGB mass loss is how the sim- 
ple molecules in the stellar photosphere grow into larger 
assemblies that form the cores for dust condensation, and 
how efficient these processes are. Little is known about the 
atmospheric chemistry, the molecular abundances within 
the dust-formation zone, and the dust-to-gas ratio in the 
wind. How do molecular chemistry and dust condensation 
depend on the stellar mass, luminosity, temperature, pul- 
sation period and amplitude, on the initial metallicity and 
on the self-enrichment of the photosphere as a result of the 
dredge-up of nucleosynthesis products? We do know that 
the transformation of some oxygen-rich AGB stars into 
carbon stars (by dredge-up on the upper AGB) changes 
the molecular and dust formation chemistry from oxygen- 
dominated to carbon-dominated, giving rise to a vastly 
different array of molecules and dust particles. Hence these 
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two types of AGB star enrich the intersteUar medium with 
completely different material. 

The Large and Small Magellanic Clouds (LMC and 
SMC) offer an excellent opportunity to study AGB mass 
loss, at metallicities of ~ 40 and 15 per cent solar, respec- 
tively. Accurate luminosities and mass-loss rates can be 
obtained in these galaxies as a result of their proximity 
and known distances of 50 and 60 kpc, respectively. Much 
work on magellanic dust-enshrouded AGB stars has been 
done (e.g., Wood et al. 1992; Zijlstra et al. 1996; van Loon 
et al. 1997, 1998). The mass-loss rate depends strongly on 
luminosity but also on temperature and pulsation prop- 
erties (Wood 1998; van Loon et al. 1999, 2005; Marshah 
et al. 2004; Whitelock et al. 2003). Dust mass-loss rates 
are lower in metal-poor environments, but this may not 
be true for the gas mass- loss rates (van Loon 2000). 

Metal-poor carbon stars are surprisingly abundant in 
C2 and C2H2, and may have a higher C/0 ratio than their 
galactic solar-metallicity equivalents (van Loon, Zijlstra 
& Groenewegen 1999; Matsuura et al. 2002, 2005). This 
may be due to the high contrast between the carbon 
self-enrichment and oxygen-poor photosphere of metal- 
poor carbon stars, and could have important consequences 
for the condensation of carbonaceous dust grains in low- 
metallicity environments. 

To investigate the formation of simple, di-atomic 
molecules such as C2 and CN, the assembly into slightly 
more complex molecules such as C2H2 and HCN, and 
the condensation into dust grains, we set out to measure 
the molecular abundances in the extended atmospheres of 
mass-losing carbon-rich AGB stars in the LMC. Because 
many of these highly-evolved objects are dust-enshrouded, 
it is impossible to perform such a study at optical wave- 
lengths. Furthermore, whilst the optical spectrum is rich 
in bands of di-atomic molecules, larger molecules are ob- 
served at infrared wavelengths. The three micron tel- 
luric window (2.9-4.1 /xm) in carbon stars shows dis- 
tinctive broad and strong absorption by C2H2 and HCN 
molecules at 3.1 /im, and often absorption by C2H2 at 
3.8 /xm and HCN at 3.57 /xm. These molecules are formed 
in the upper atmospheric layers by hydration of the di- 
atomic molecules that are seen at optical wavelengths. 
C2H2 is a possible building block for Polycyclic Aromatic 
Hydrocarbons (PAHs) (Latter 1991), and must play a key 
role in the condensation of carbonaceous grains in H-rich 
environments (Gail & Sedlmayr 1988). The HCN band 
strength may be a measure of the nitrogen abundance. 

Here wc present ESO/VLT three micron spectra of a 
large sample of heavily dust-enshrouded carbon stars in 
the LMC, and an analysis of their molecular atmospheres. 

2. Observations 

2.1. Targets 

The targets for our 3-4 /xm spectroscopy were drawn from 
a sample of mid-IR sources detected with the Mid-course 
Space experiment (MSX) (Egan, Van Dyk & Price 2001; 



Wood & Cohen 2001). Our principal objective was to 
establish the nature and chemical composition of heav- 
ily dust-enshrouded stars for which optical spectroscopy 
would not be possible, to help select targets for our Spitzer 
Space Telescope programme #3505 (PL P.R. Wood). 
Hence wc selected objects by excluding all objects that 
have optical photometry in the MACHO database, bright 
objects with a Kg-band magnitude < 9 mag, and objects 
that are not extremely red with a J-Kg colour < 5 mag. 
The actually observed objects are mostly of 8"^ or 9**^ 
magnitude in the 3-4 jim range. Throughout the paper 
we refer to the MSX-selected sample by their MSX LMC 
catalogue numbers. 

Six more dusty objects were observed that were known 
to be in or near the populous LMC clusters NGC 1903, 
NGC 1984 and NGC 1978 (van Loon, Marshah & Zijlstra 
2005). VLT/ISAAC spectra covering the same spectral re- 
gion are available from the literature for nine more IR 
objects in the LMC (Matsuura et al. 2005). Of these, 
IRAS 04496-6958, IRAS 04557-6753 and LI-LMC1813 
are cluster members (van Loon et al. 2005b). The IR pho- 
tometric data of the cluster objects can be found in van 
Loon et al. (2005b) and are not repeated here. 

2.2. Spectroscopy 

The ISAAC instrument on the European Southern 

Observatory (ESO) Very Large Telescope (VLT) "Antu" 
at Paranal, Chile, was used on the nights of 6, 7 and 8 
December 2003 to obtain low-resolution long-slit spectra 
between 2.9 and 4.1 /im. At a slit width of 2" the spectral 
resolving power was essentially determined by the seeing 
and was 300^i? ^ 700. The smallest attained resolution 
element is sampled by three pixels on the Aladdin array. 

The spectra were obtained by chopping and nodding 
with a throw of 10" (15" in some cases to avoid chopping 
onto another star) to cancel the high background at these 
thermal-IR wavelengths, and jittering within 2" to remove 
the effects of bad pixels. This produced three spectra on 
the final combined frame, with the central spectrum twice 
as bright as the other (inverse) spectra. Exposure times 
varied between 12 and 60 minutes. After removing cos- 
mic ray hits and correcting the spectra for their small 
inclination with respect to the array, all spectra were ex- 
tracted using an optimal extraction algorithm. The ex- 
tracted spectra were then aligned and combined. Spectra 
of an internal Xe-|-Ar lamp were used to calibrate the 
wavelengths; the spectra are sampled on a grid of 14 A 
elements. 

The relative spectral response calibration was obtained 
by dividing by the spectrum of the B-type standard 
star HIP 020020. This removed most of the telluric ab- 
sorption lines — except the saturated methane line at 
3.32 /xm — but introduced artificial emission features 
due to photospheric absorption lines in the spectrum of 
HIP 020020. Each spectrum was therefore multiplied by a 
synthetic spectrum appropriate for HIP 020020, in which 
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Table 1. List of targets, in order of increasing Right Ascension (2MASS coordinates, in J2000). MSXLMC (Egan, van 
Dyk & Price 2001; targets for Spitzer Space Telescope programmes are in boldface) and IRASPSC designations are 
given where available, together with ncar-IR magnitudes (1=2MASS, 2=MSSS0, 3=spcctroscopy acquisition images) 
and MSX band- A and IRAS 12 and 25 /xm flux densities (in mJy). Values marked with a colon are suspect. References 
to the identification as a dust-enshrouded AGB star are as follows: l=Egan, Van Dyk & Price (2001); 2=Loup et 
al. (1997); 3=Reid, Tinney & Mould (1990); 4=Trams et al. (1999); 5=van Loon et al. (1997); 6=van Loon (2000); 
7=Zijlstra et al. (1996); 8=this work. 
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LMC oxygen-rich stars from Matsuura et al. (2005) 
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Spitzer programmes: a = GTO #200 (Houck); b = DDT #1094 (Kemper); c = GO #3426 (Kastner); d = GO #3505 (Wood). 
Notes: A = SAAO monitoring at 3.5 /im (Whitelock et al. 2003); B = BMB-R46, SHV 0521050-690415; C = HV916. 



wc have inserted absorption lines due to Bra 4.052, Pf7 
3.741, PfJ 3.297, Pfe 3.039 and PfC 2.873 iim of depths 
sufScient to remove the spurious emission lines. Small 
wavelength shifts and differences in effective resolution 
can cause imperfections in the correction. Fig. 1 shows 
the observed uncorrected spectra of HIP 020020 and tar- 
get star MSX LMC 47, a low-resolution Kurucz spectrum 
(Teff = 18,000 K, log 5 = 4.0 and [Fe/H]=0) and our syn- 



thetic spectrum. The observed spectra show the extent to 
which the telluric absorption affects the spectrum espe- 
cially around 3.2-3.3 ^m, and the effective band width 
between 2.85 and 4.15 /xm. 

A K-a clipping algorithm was applied to remove spikes 
from the spectrum that are caused by hot pixels. Care 
was taken to ensure that astrophysical features were not 
affected by this procedure. 
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X [miti] 

Fig. 1. Uncorrected ESO-VLT/ISAAC spectra of the 
B-type standard star HIP 020020 and a target star, 
MSX LMC 47, plus a low-resolution Kurucz spectrum and 
our synthetic spectrum appropriate for the photosphere of 
HIP 020020. 

The conditions were very similar for all three nights. 
The seeing as measured on the acquisition images was 
around 0.5" most of the time but occasionally worsened 
up to about 1" (e.g., when observing at high airmass of 
~2). The relative humidity and air temperature were fairly 
constant at 8-18 per cent and 12-15°C, respectively. 

2.3. Photometry 

IR photometry (Table 1) is compiled for all targets and 
for the Matsuura et al. (2005) objects, except for cluster 
members (see van Loon et al. 2005b). A few have been 
monitored at near-IR wavelengths to measure the pulsa- 
tion period and amplitude (Table 2). 

JKL-band photometry for the targets was obtained at 
Mount Stromlo and Siding Springs Observatory (MSSSO) 
as part of a long-term monitoring campaign (Wood et al., 
in preparation). In addition, Kg-band photometry from 
the 2-Micron All-Sky Survey (2MASS) is available for al- 
most all targets, but most targets are too faint in the 
J-band and many are even too faint in the H-band to 
have 2MASS data at those wavelengths. Very few targets 
are detected in the DENIS survey, which is inferior to 
2MASS in terms of sensitivity. Only three MSX-selected 
targets have a DENIS Kg-band measurement (Tabic 2). 
Additional L'-band photometry is available from the spec- 
troscopy acquisition images, accurate to within a few 0.1 
mag. 

All prime targets have reliable MSX band A (8.3 fim) 
flux densities listed in the MSX Point Source Catalogue 



(Version 2.3). We also collected InfraRed Astronomical 
Satellite (IRAS) scans from the IRAS data server^ and 
used the Groningen Image Processing SYstcm (GIPSY) 
software with the SCANAID tool to reconstruct a cut 
through the emission on the exact position. This is supe- 
rior to the IRAS Point Source Catalogue values and works 
well for isolated sources for which often reliable flux den- 
sities can be estimated down to a level of a few 0.01 Jy 
at 12 and 25 fim. A few targets have also been observed 
with the Infrared Space Observatory (ISO) (Table 2). 

3. Results 

3.1. IR carbon stars 

The majority (28) of our targets turn out to be carbon 
stars, on the basis of strong carbonaceous molecular ab- 
sorption features seen in the spectrum around 3.1 and of- 
ten 3.8 /xm (Fig. 2). Out of these, 10 are part of Spitzer 
Space Telescope GO programme #3505 (PI: P. Wood). 

It is worth noting that MSX LMC 307 has the longest 
known pulsation period of any carbon star: 889 d (Wood 
1998) or 939 d (Whitclock et al. 2003). 

MSX LMC 55 is misidentified in the literature with the 
M-type giant HV 11977, of which we also took a spectrum. 

Another mistaken identity occurs for MSX LMC 45: 
Leisy et al. (1997) identified the IRAS source with the low- 
excitation Planetary Nebula SMPLMC36 (Sanduleak, 
MacConnell & Philip 1978), located at 9" WNW from 
an extremely red 2MASS source. Our 3-4 /xm spectrum is 
of the 2MASS object, whilst the PN is a target for Spitzer 
Space Telescope GTO programme #103 (P.I.: J. Houck). 
Although the K^-L' colour of SMP LMC 36 (Table 3) is 
very red at 2.76 mag, MSX LMC 45 is much redder at 4.95 
mag and with L' = 7.87 mag much brighter in the thermal 
IR than SMP LMC 36 which has L' = 11.24 mag. For ah 
eff'ects the emission detected at 8.3 (MSX), 12 and 25 /xm 
(IRAS) arises from the carbon star. 

3.2. IR objects of unknown type 

For two of our targets, the 3-4 /xm spectrum does not lead 
to an unambiguous photospheric chemical classification 
(Fig. 3). Coincidentally, these arc the only two objects in 
our target sample that are part of Spitzer Space Telescope 
programmes other than GO #3505. Three more objects 
turned out to be massive Young Stellar Objects; these are 
analysed in detail in van Loon et al. (2005c) and Oliveira 
et al. (in preparation). 

The 3-4 fim spectrum of MSX LMC 196 displays a hint 
of absorption around 3.0-3.1 and 3.7-3.8 /xm, in which 
case it is a carbon star. The molecular features may be 
heavily veiled by the strong dust continuum emission. 
Indeed, with F12 = F25 its [12]-[25] colour is redder than 
any of the IR carbon stars. 

On the basis of the IR colours, van Loon et al. (1998) 
classified the dust around MSX LMC 733 as oxygen-rich, 

^ htt p : / / www. astro .rug.nl /IRAS- Server / 
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Table 2. Additional data from DENIS and ISO, and pulsation properties from near-IR monitoring. The flux densities 
are in mJy, pulsation periods in days and Kg-band amplitudes in magnitudes. The references are as follows: l=Loup 
et al. (in preparation); 2=Trams et al. (1999); 3=van Loon et al. (2005a); 4=van Loon et al. (2005b); 5=Whitelock et 

al. (2003): 6=Wood (1998). 
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67 




51 






541 


0.65 


1,2,5 


LMC oxygen-rich stars from 


Matsuura et al. (2005) 


















IRAS 05128-6455 


17.87 


13.43 






226 


61 


PHOT-S 


409 


0.30 


2,3,5 


IRAS 05148-6730 


10.47 


7.34 






380 


176 




951 


0.24 


2,3,5 



Table 3. Serendipitous objects, listed in order of increasing Right Ascension (2MASS coordinates, in J2000), together 
with near-IR magnitudes from 2MASS (JHKg) and the spectroscopy acquisition images (L'). The PN candidate 
and Spitzer Space Telescope target SMPLMC36 is included as we obtained L'-band photometry for it, although no 
spectrum was taken. 



Name 




5 (° ' ") 


J 


H 


K, 


L' 


Type 


MSXLMC 50 B 


5 02 16.9 


-68 05 24 


11.25 


10.40 


10.17 


9.87 


? 


HV 11977 


5 05 45.9 


-68 38 54 


10.57 


9.77 


9.53 


9.24 


M giant 


SMP LMC 36" 


5 10 39.7 


-68 36 05 


16.25 


15.61 


14.00 


11.24 


Planetary Nebula? 


MSXLMC 349 B 


5 17 24.4 


-68 55 06 


11.62 


10.50 


9.83 


9.98 


Carbon star 


HV2532 


5 26 27.4 


-69 10 56 


8.93 


8.08 


7.71 


7.42 


M4 supergiant 


MSXLMC 1780 B 


5 32 59.0 


-68 23 51 


11.10 


10.54 


10.46 


10.12 


Early-type star? 



Spitzer programme: a = GTO #103 (P.I.: J. Houck). 



but in a re-analysis including ISO data Trams et al. (1999) 
classified it as carbon-rich. They also present a CAM- 
CVF spectrum, in which broad emission around 11.3 /xm 
was interpreted as SiC, which would confirm its carbon- 
rich nature. However, its CAM-CVF spc;c;truni is not un- 
like that of OH/IR stars where the silicate feature is 
seen in self-absorption. Overall the 3-4 ^m spectrum of 
MSXLMC 733 is similar to the ISO-SWS spectrum of the 
M8 giant CDGru (Vandenbussche et al. 2002). 



3.3. Serendipitous objects 

In several cases a spectrum could be taken of another IR- 
bright object, simultaneously with the principal target. 
These "serendipitous" objects (Fig. 4 & Table 3) are not 
generally dust-enshrouded. 



The M4 supergiant HV2532 is responsible for most 
of the mid-IR emission associated with IRAS 05267—6913 
{Fi2 ~ 0.20 Jy, F25 ~ 0.25 Jy). The MSX data clearly 
separates its emission, MSX LMC 592 [Fa = 230 mJy), 
from that of our spectroscopy target, MSX LMC 774. 

3.4. Cluster IR objects 

We took spectra of six cluster IR objects from van Loon et 
al. (2005b) (Fig. 5). Four of these, NGC 1903-IRl and IR2 
and NGC 1978-IRl and IR2 are carbon stars. The other 
two, NGC 1984-IRl and IR2 arc oxygen-rich objects; IRl 
(=IRAS 05280-6910) is a source of OH (Wood et al. 1992; 
Marshall et al. 2004) and H2O maser emission (van Loon 
et al. 2001) whilst IR2 is the Ml supergiant WOH G347. 

Two objects in our MSX-selected target list, 
MSX LMC 484 and 654, and three of the carbon stars 
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3.5 



Fig. 3. ESO/VLT 3-4 /xm spectra of IR objects of un- 
known type. MSXLMC196 shows a hint of absorption 
around 3.1 /xm and 3.65-3.80 fim, and may thus be a heav- 
ily veiled carbon star. The overall shape of the spectrum 
of MSX LMC 733 resembles that of cool oxygen-rich stars 
better. 



with 3-4 iim spectra published in Matsuura et al. (2005), 
IRAS 04496-6958, IRAS 04557-6753 and LI-LMC1813 
were also found to be associated with clusters (van Loon 
et al. 2003; van Loon et al. 2005b). 

4. Analysis 

4.1. Luminosities and mass-loss rates 

As a result of the selection criteria, the MSX-selected tar- 
gets have very red Kg-L' (Fig. 6). This is caused by cir- 
cumstellar extinction diminishing the Ks-band brightness 
and circumstellar dust emission enhancing the L'-band 
brightness. Hence their faint Kg-band brightness is not a 
direct consequence of their bolometric luminosity. IR car- 
bon stars are objects which are known to be variable with 
amplitudes up to 2 magnitudes at near-IR wavelengths 
(e.g., Wood 1998; Whitelock et al. 2003). Although the 
2MASS, MSX and IRAS data refer to combined measure- 
ments taken at different epochs, the groundbased near-IR 
data that we use here are single-epoch measurements. 

The spectral energy distributions of the IR carbon 
stars were modelled with the radiative transfer code 
DUSTY (Ivezic, Nenkova & Elitzur 1999). The density dis- 
tribution is based upon a hydrodynamic computation of 
a dust-driven wind at constant mass-loss rate. The model 
was then scaled to match the overall observed Spectral 
Energy Distribution (SED), knowing the distance to the 
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0.5 

2 

1.5 
1 
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H 1 1 1- 



H 1 1 1- 



MSX LMC 349 B 
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MSX LMC 1780 B 




3 3.5 X [fj,m] 4 

Fig. 4. Serendipitous ESO/VLT 3-4 ^m spectra (see 
Table 3) . They are all oxygen- rich except MSX LMC 349 B 
which is an optically bright (yet anonymous) carbon star. 



LMC (50 kpc), which then yields an accurate measure- 
ment of the bolometric luminosity. 

A blackbody was used to represent the underlying stel- 
lar radiation field, with a temperature of Tcs = 2500 
K. We used 0.1 /xm grains of amorphous carbon dust 
(Henning & Mutschke 1997). To obtain the mass-loss rate 
one has to know the dust grain density, for which wc adopt 
Pgrain = 3 g cm~^, and the gas-to-dust mass ratio, for 
which we adopt Pgas/pdust = 500 at the metallicity of the 
LMC (two to three times below solar). The results are 
summarised in Table 4. 

We usually found that a satisfactory match between 
the model and the observed SED required a rather low 
temperature for the dust at the inner edge of the envelope, 
Tdust ~ 600 to 700 K. The derived bolometric luminosity 
and mass-loss rate do not however depend sensitively on 
the inner-edge dust temperature. The estimated internal 
accuracies of the luminosity and mass-loss rate are about 
5 and 10 per cent, respectively. 

The mass-loss rates axe correlated with the bolometric 
luminosities (Fig. 7). Most stars straddle the classical limit 
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Fig. 2. ESO/VLT 3-4 /rni spectra of MSX-selected IR carbon stars. The spectra have been normaUsed to i^s.s = 1. 



to the mass-loss rate for single scattering (long-dashed line 
in Fig. 7; van Loon et al. 1999b). The observed sequence 
is steeper than the classical limit though, with luminous 
stars experiencing mass-loss rates that are higher than this 
limit — but not as high as the maximum mass-loss rates 
observed in the LMC (short-dashed lines in Fig. 7; van 
Loon ct al. 1999b). The two carbon stars with particu- 
larly long pulsation periods, P > 800 d, and amplitudes, 
AKs > 1.7 mag, MSX LMC 307 and IRAS 05112-6755 
(Table 2) are amongst the stars with the highest mass- 
loss rates and luminosities. 



On the whole the mass- loss rates of the IR carbon stars 
are higher than those of optically bright carbon stars (van 
Loon et al. 2005a), but amongst the IR carbon stars them- 
selves are stars that have relatively low mass-loss rates for 
their high luminosity. Four of the carbon stars from van 
Loon et al. (2005b) and Matsuura et al. (2005) do not ap- 
pear in Fig. 7 as they have significantly lower mass-loss 
rates than the stars in the MSX-selected sample. These 
stars may bc^ in a slightly different evolutionary phase, 
for instance if they have not yet entered the superwind 
stage or if they pass through a more quiet episode in their 
thermal-pulsing cycle. Differences in mass may also con- 
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Fig. 2. (Continued) 



tribute to a spread in the mass-loss rate-luminosity plane 
— see van Loon et al. (2005b) for a discussion of the 
mass-loss rates and luminosities of cluster IR objects in 
the Magellanic Clouds. The by far most luminous carbon 
star that we analyse here, IRAS 04496—6958 also lies be- 
yond the boundaries of the diagram. Its mass-loss rate is 
fairly modest for its high luminosity. In this case there are 
strong indications for a relatively massive progenitor star 
(Afinitiai ~ 5 Mq) that may have experienced a thermal 
pulse recently (van Loon et al. 2005b). For the purpose of 
the spectral analysis, we divide the sample of IR carbon 



stars into four groups on the basis of their mass-loss rate 
and luminosity (Fig. 7 & Table 4). 



Stars with high mass-loss rates have a red Kg-L' colour 
or, if having a relatively modestly red Kg-L' colour they 
have a red L'-[8.3] colour (Fig. 8). The fact that the L'- 
[8.3] colour shows some correlation with the Kg-L' colour 
but with a large scatter implies that besides the optical 
depth of the dust shell other parameters play a role. These 
may include the dust temperature at the inner edge of 
the envelope and perhaps the dust type, or photometric 
variability. 
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Fig. 5. ESO/VLT 3-4 nm spectra of additional objects in or near populous star clusters (van Loon et al. 2005b). They 
are all carbon stars except NGC 1984-IRl and -IR2 which are an OH/IR star and an optically bright red supergiant, 
respectively. 
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Fig. 6. Absolute Kg brightness versus Kg-L' colour- 
magnitude diagram of the MSX sample, where the ob- 
jects for which we present spectra (Table 1) are indicated 
by squares and circles. 

4.2. Identification of spectral features 

A grand average of all spectra of the IR carbon stars (Fig. 
9) shows all the spectral features of interest at a much 
reduced noise level. Examples of four types of spectral 
appearance are given in Fig. 10, to show the similarities 
and differences amongst the IR carbon stars. 
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Fig. 7. Mass-loss rate versus luminosity diagram for IR 
carbon stars, which we use to subdivide the sample into 
four groups represented by different symbols. Stars from 
the cluster sample of van Loon et al. (2005b) and previ- 
ously published spectra of Matsuura et al. (2005) are kept 
separate (crosses). 



Hydrogen absorption lines of the Brackett and Pfund 
series are seen in many of the spectra, but this may be 
due to residuals in the correction for the same lines in the 
telluric standard star spectrum. There is no evidence in 
any of the carbon stars for strong hydrogen line emission 
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Table 4. Bolometric luminosities, mass-loss rates, inner- 
edge dust temperatures and optical depth of the IR carbon 

stars. They arc grouped (last column) as follows: l=high 
M and L; 2=moderate M and L; 3=low M and L; 4=low 
M but high L. 



Name 


1 /T N 

log(L) 
(L©) 


log(M) 
{Mq yr" 


Tdust 

) (K) 


TV 


G 


IR carbon stars 












MSX LMC 1042 


3.91 


—4.80 


600 


13 


3 


MSXLMC 1249 


4.07 


-4.84 


700 


10 


4 


MSXLMC 50 


3.95 


—4.51 


650 


33 


2 


MSX LMC 91 


3.85 


—4.77 


650 


19 


3 


MSXLMC 93 


3.88 


—4.82 


700 


17 


3 


TV JTl^^T" X TV JT 

MSX LMC 55 


4.04 


—4.74 


450 


4 


4 


MSXLMC 45 


4.31 


-4.30 


650 
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MSXLMC 47 


4.07 


—4.51 


650 


25 


2 


TV X TV IT C\-t r\ 

MSXLMC 219 


4.03 


—4.58 


700 


25 


2 


TV !rC~^'\T X TV T ■ "x .1 ■ -\ 

MSX LMC 232 


3.90 


—4.73 


600 


17 


3 


MSX LMC 223 


4.00 


—4.54 


700 


30 


2 


TV Tn"\7" X TV T ir\r\r\ 

MSXLMC 202 


3.78 


—4.76 


700 


25 


3 


MSXLMC 225 


4.01 


—4.48 


650 


32 


2 


MSXLMC 221 


3.94 


—4.65 


550 


17 


2 


TV TrT\7" X TV T O A r\ 

MSXLMC 349 


3.90 


—4.54 


700 


37 


2 


MSXLMC 307 


4.15 


-4.38 


650 


31 


1 


TVTn"t7"X TVT i\ A t 

MSXLMC 341 


3.80 


—4.75 


700 


25 


3 


A/rGV T A/T/^ AQA 

MbA LiiViU 484 


6.1 L 


—4. (0 


iUUU 


oi 


Q 
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MSX LMC 774 


4.00 


-4.60 


600 


20 


2 


MSX LMC 677 


3.88 


-4.85 


700 


15 


3 


MSXLMC 635 


4.27 


-4.31 


650 


29 


1 


MSXLMC 692 


4.01 


-4.70 


700 


18 


4 


MSXLMC 644 


4.02 


-4.71 


700 


17 


A 


MSXLMC 654 


3.88 


-4.67 


700 


20 


o 
z 


MSXLMC 1780 


4.07 


-4.35 


600 


37 


1 


MSX LMC 743 


3.76 


-4.85 


700 


20 


6 
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3.94 


-4.61 


700 


28 


o 
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4.07 


-4.69 


700 


16 


4 


Cluster carbon stars from van Loon ( 


t al. (2005b) 
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-5.11 
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3.84 


-5.43 
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3.73 


-5.05 
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-4.54 


720 
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700 
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3.99 


-4.57 
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IRAS 05112-6755 


4.22 


-4.25 
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BMB-R46 


4.19 


-5.27 


700 


1 




LI-LMC 1813 


4.09 


-4.33 


600 


28 





arising from shocks in their pulsating atmospheres. A few 
spectra show an array artifact causing a sharp feature in 
the extracted spectrum around 3.99 /xm. 

The strongest absorption band, between about 2.95 
and 3.15 fim, is seen in all IR carbon stars, and is due 
to a blend of HCN and C2H2. The shape is not uni- 
versal, although often with characteristic peaks at 3.01 
and 3.11 /^m it sometimes has a rather triangular shape 
peaking at 3.01 /xm only (e.g., MSXLMC 1249). Often, 
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Fig. 8. L'-[8.3] versus Kg - L' colour-colour diagram of the 
IR carbon stars. Redder colours are generally indicative 
of a larger optical depth of the circumstellar dust enve- 
lope, but there is significant scatter due to other (cir- 
cum)stellar parameters and non-contemporaneous pho- 
tometry. Symbols are as in Fig. 7. 

the absorption seems to extend to as far as 3.5 /im (e.g., 
MSXLMC 221), but it is not certain that this is still part 
of the same band. Alternatively, the HNC isomer of HCN 
is expected to contribute in precisely this spectral region 
(Harris et al. 2003). 

Most IR carbon stars exhibit a broad depression in 
the spectrum around 3.8 /im. In the LMC, the absorption 
almost always peaks around 3.75 /im, whereas in galac- 
tic carbon stars absorption is normally seen around 3.9 
/im instead. In our sample, only the less extreme object 
NGC 1903-IRl shows evidence for distinct absorption at 
3.9 /im. CH was suggested to explain a series of absorption 
lines around 3.75 /im in the spectra of the galactic carbon 
stars TXPsc and WZCas (Aoki, Tsuji & Ohnaka 1998), 
but the CH does not lead to a broad absorption band. 
Instead, the identification with C2H2 by Cernicharo et al. 
(1999) matches the broad absorption seen around 3.75 /im 
in the galactic dust-enshrouded carbon star IRC -1-10216. 
HCN is held responsible for the 3.9 /im component (Aoki 
et al. 1998; Harris et al. 2003), but CS may also contribute 
(Aoki et al. 1998). 

A few quite narrow absorption features are seen, su- 
perimposed onto the 3.8 /im absorption band (Fig. 9). The 
3.70 /im feature is the most conspicuous of these, but it is 
always seen in conjunction with a similar feature at 3.76 
/im. A weaker feature at 3.78 /im may be related, and 
there seems to be another feature that is however blended 
with the Pf7 line at 3.74 /im. Aoki et al. (1998) interpret 
3.70 and 3.77 /im features as being part of the CH(3-2) 



van Loon et al.: 3 jjm spectra of dust-enshrouded stars in the LMC 



11 




3 3-5 A [^m] 4 

Fig. 10. Examples of four types of spectral appearance 
amongst the IR carbon stars. 

fundamental band, but the stronger associated 3.64 and 
3.67 /xm features (cf. J0rgensen et al. 1996) are not seen in 
our spectra (whilst they are seen in TX Psc and WZ Cas) 
— except perhaps in the enigmatic object MSXLMC 196 
(Fig. 3). The features might be part of the C2H2 band 
instead (Cernicharo et al. 1999). 

CH(l-O) may be responsible for the series of narrow 
absorption lines seen around 3.3-3.5 fim in the spectrum 
of MSXLMC 349 B (Fig. 4). The spectrum of this optically 
bright carbon star is in fact slightly depressed up to a 
wavelength of about 3.75 /im. 

A comparison between the average spectra for the four 
groups in terms of mass-loss rate and luminosity (Fig. 
11) fails to reveal any striking correlation of the spectral 
features with either of these parameters. The dispersion 
within each group average is quite large, but the error in 
the mean difference between the group average and the 
grand average indicates statistically significant differences 
between the groups. The most obvious difference is in the 
slope of the spectral continuum, which depends on the 
optical depth of the circumstellar dust envelope: the con- 
tinuum is redder at higher mass-loss rates, but at the same 
mass-loss rate it is redder at lower luminosity. 

4.3. Molecular band strengths 

To investigate the band strengths of some of the spectral 
features, we compute equivalent widths (Tables 5 & 6). 
For the narrow 3.57 and 3.70 fini bands it is important to 
realise that the band definitions are appropriate for our 
spectra and hence include the velocity shift of the LMC 
(0.003-0.004 fim) and the internal wavelength calibration 
accuracy (±0.002 fim) . In other data sets these definitions 




Fig. 11. Average spectra for the four groups of stars se- 
lected on the basis of their mass-loss rate and luminosity, 
plotted together with their difference spectra with respect 
to the grand average of all IR carbon star spectra (see Fig. 
9). The dots indicate the dispersion in the group averages 
and the error in the mean difference between the group 
and grand averages. 



Table 5. The band and blue/red continuum definitions for 
the computation of equivalent widths. Wavelengths are in 
/Ltm. 



Feature 


band 


blue continuum 


red continuum 


3.1 


2.950-3.350 


2.900-2.950 


3.350-3.400 


3.57 


3.560-3.580 


3.510 3.560 


3.580-3.630 


3.70 


3.700-3.715 


3.680-3.700 


3.715-3.720 


3.8 


3.600-4.000 


3.500-3.600 


4.000-4.100 



might therefore need a slight readjustment. We recom- 
puted the equivalent widths at 3.1, 3.57 and 3.8 fj,m in the 
spectra from Matsuura et al. (2005) after converting their 
F^ units to our Fa units. Different units cause different 
spectral slopes and hence different continuum definitions, 
which affects the equivalent width measurements. We also 
made a slight adjustment to the definition of the blue edge 
of the 3.1 yum band. 

The uncertainties in the equivalent width values are 
obtained from the standard deviations of the values for 
the eight oxygen-rich objects: a = 170, 2.23, 3.03 and 
175 A in the equivalent widths at 3.1, 3.57, 3.70 and 3.8 
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Fig. 9. Grand average of all IR carbon star spectra (Table 1), with identifications of the main features discussed in 
the text. 



/zm, respectively. This assumes that oxygen-rich objects 
have zero equivalent widths, which is not strictly true; the 
quoted uncertainties are therefore conservative estimates. 
The 3.57 and especially the 3.70 /zm features are quite 
difficult to measure in individual spectra. 

All four groups of IR carbon stars show a very clear 
correlation between the strength of the 3.8 ^m band and 
that of the 3.1 fim band (Fig. 12), with a Pearson correla- 
tion coefficient of p = 0.80: an increase in column density 
of C2H2 would make both bands appear stronger. The 3.1 
yum band is always prominent, within a range in equivalent 
width of a factor two. The 3.8 /im band is always seen in 
conjunction with strong 3.1 fini absorption, W3.1 > 1000 
A {p = 0.85, including the additional stars — crosses in 
Fig. 12), but it is not always visible in spectra with weaker 
3.1 fim absorption, when the 3.1 /zm band may be domi- 
nated by the contribution from HCN. 

For W3.8 < 450 A; the correlation between the 3.70 and 
3.8 /um band strengths is much better {p = 0.41) than that 
between the 3.57 and 3.8 /xm band strengths {p = 0.11). 
This might imply that the 3.70 ^m band is associated with 
C2H2. At W3.8 > 450 A veiling and changes in excitation 
conditions (see below) may become important and affect 
the correlations. 

Some of the scatter in the equivalent width distribu- 
tions is intrinsic. For instance, MSXLMC91 exhibits ex- 
tremely strong 3.57 and 3.70 /im bands but only moder- 
ately strong 3.1 and 3.8 iim bands, whilst MSXLMC221 
shows exactly the opposite behaviour. Another extreme 
object is the relatively blue MSXLMC 1249, whoso 3.1 ^m 
band resembles that of water ice (which is not expected to 
exist in a warm carbon-dominated chemistry) and which 
does not exhibit any 3.70 or 3.8 /zm absorption although 
it does show a small dip at 3.57 /zm. The differences be- 
tween extreme examples of 3-4 /xm spectra (Fig. 10) are 



much greater than the subtle differences as a function of 
luminosity and mass- loss rate (Fig. 11), suggesting that 
other parameters influence the shape and strength of the 
molecular bands. 

The shape of the spectrum between about 3.1 and 3.7 
/j.m can be dramatically different for different carbon stars 
(see Fig. 10), without it being clear what absorption band 
— if any — is responsible for this. To investigate this 
further we compute the spectral slopes, s(A), around 3.15, 
3.24, 3.43 and 3.57 /zm within the regions 3.12-3.18, 3.20- 
3.28, 3.36-3.50 and 3.50-3.64 (excluding 3.56-3.58) /zm, 
respectively. We then compute the pair-wise differential 
spectral slopes, s'(A) (Table 7), to quantify the degree 
to which the spectrum displays a "knee" such as is seen, 
for example, at 3.2 /zm in MSX LMC 307 or at 3.5 /zm in 
MSXLMC 221. If there is no change in spectral slope the 
differential spectral slope is zero. 

The differential slopes at 3.2 and 3.5 /zm show a strong 
anti-correlation (Fig. 13; p = —0.57, ov p = —0.62 if the 
noisier spectra of group 3 sources are omitted). There is a 
gradual transition between two extreme types of spectral 
shape: those with a narrow 3.1 /zm band and a flat con- 
tinuum between 3.2 and 3.7 /zm ("cup"), and those with 
absorption extending from the 3.1 /tm band all the way 
until 3.5 /zm where it meets the continuum or dips into the 
3.8 /zm band ("bowl"). Interestingly, the best examples of 
the "cup" type are all moderate to high mass-loss rate 
and luminosity objects (including IRAS 04286— 6937) al- 
though the "bowl" type also occurs amongst these classes 
of stars. 

The strength of the 3.70 /tm feature is correlated with 
the differential slope at 3.2 /zm (Fig. 14; p = 0.45). This 
suggests that whatever causes the broad absorption be- 
tween 3.2 and 3.5 /zm is anti-correlated with the 3.70 /tm 
feature. On the other hand the strength of the 3.8 /tm 
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band is correlated with the differential slope at 3.5 /xm 
(p = 0.72) and mildly anti-cor related with that at 3.2 jim 
(p = —0.33), suggesting a link between the 3.2 -3.5 and 
3.8 /Ltm absorptions. We interpret this as an effect of dif- 
ferent excitation conditions. A bowl-shaped 3.1 /um band 
indicates a high excitation temperature, which would then 
also result in broad absorption by C2H2 at 3.8 i^m. A low 
excitation temperature would result in narrow absorption 
features of C2H2, between 3.70-3.78 /um. 



4.4. Relation between the molecular atmosphere and 
circumstellar dust envelope 

Can some of the differences in the molecular bands be 
reflected in differences in the circumstellar dust envelope? 
Stars with the highest mass-loss rates and luminosities 
always have strong 3.8 /xm bands and almost always do 
they also have a conspicuous 3.70 fxm feature (Fig. 15). 
Stars with low mass-loss rates of M < 10~^ Mq yr~^ 
tend to have weak or absent 3.70 and 3.8 /xm absorptions 
whilst the 3.57 /um band is often visible in these stars. 

At high mass- loss rates of M > 2.5 x 10"^ M© yr'^ the 
3.1, 3.57 and 3.8 /xm absorption bands increase in strength 
with increasing mass-loss rate. However, at lower mass- 
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Table 6. Equivalent widths (in A). The chemical classi- 
fication into carbon stars (C) and oxygen-rich stars (o) is 

based on the presence or absence of the 3.1 /im feature, 
supplemented by information from the literature. 
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Table 7. Differential spectral slopes s'(3.2) = s(3.15) — 
s(3.24) and s'(3.5) = s(3.43) - s(3.57) (see text), and the 
fractional contributions of the dust emission, for carbon 
stars. 
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loss rates the bands are generally stronger than around 
3 X 10~^ M© yr~^, albeit with considerable spread. The 
same trends are seen for the 3.70 /zm band. This prob- 
ably reflects a competition between the effect of veiling 
(reducing the equivalent width) and an optically thicker 
atmosphere in connection with a higher mass-loss rate (in- 
creasing the equivalent width). 

Dust column densities increase with increasing mass- 
loss rate but also with decreasing luminosity as low- 
luminosity objects are more compact. Molecular column 
densities may be affected in a similar way, and the dust 
optical depth might therefore show a stronger correlation 
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Fig. 13. Change in slope around 3.5 /xm versus that 
around 3.2 /zm; there is a clear correlation. Symbols are 
as in Fig. 7. 
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Fig. 14. Change in slope around 3.2 and 3.5 /ttm versus 
the equivalent width of the 3.70 and 3.8 /im bands. The 
3.70 /xm band is generally stronger for larger S3 2 and 
smaller S3 5, whilst the 3.80 /xm band shows the opposite 
behaviour. Symbols are as in Fig. 7. 



with the molecular column density than either the mass- 
loss rate or luminosity would do. The 3.1 and 3.8 /xm bands 
become weaker at progressively larger dust optical depth 
(Fig. 16). This is likely due to veiling by the increasing 
contribution of dust emission. The decrease in equivalent 
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Fig. 15. Equivalent widths versus mass- loss rate. There 
is no clear correlation other than that the absorption is 
generally strong at the highest mass-loss rates. Symbols 
are as in Fig. 7. 



width is indeed larger at 3.8 /xm than at 3.1 /itm. The 3.57 
and 3.70 /im equivalent width seems to remain constant, 
though, which may suggest that they would be seen to 
increase with optical depth if the effects of veiling were to 
be removed. At relatively modest optical depth, tv < 10, 
the 3.57 iJ,m band is often quite strong but the 3.70 /xm 
band is often absent. 

The effects of veiling by dust emission can in principle 
be removed, as dusty returns the fractional contribution 
of the dust emission to the total flux density as a function 
of wavelength. We tabulate these values at wavelengths of 
3.0, 3.5 and 4.0 /xm (Table 7). Unfortunately, these frac- 
tions are ill-determined: they depend sensitively on the 
amount and temperature of the warmest dust, which is 
also the most time- variable component. Hence correcting 
the equivalent widths in this way increases the overall scat- 
ter in a sub-set of the sample of carbon stars (Fig. 17), 
which is especially due to stars with fairly modest mass- 
loss rates. However, for at least half the sample of stars the 
scatter is reduced and the 3.1 /xm equivalent widths of the 
four groups of stars are more in line with each other (Fig. 
17). This suggests that veiling is the main driver behind 
variations of the 3.1 /um equivalent width with mass-loss 
rate and luminosity. 

Consequently the photospheric 3.1 /zm equivalent 
width seems to depend very little on stellar or circum- 
stellar parameters, and we propose that this band may 
be saturated for the MSX-selected sample of IR carbon 
stars. Comparing the dust continuum contribution (Table 
7) with the maximum depth of the 3.1 fim absorption in 
the spectra of Fig. 2, it appears that part of the absorp- 
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Fig. 16. Equivalent widths versus dust optical depth. The 
error bars refer to the mean values and standard devia- 
tions in the intervals ry G (0,10], (10,20], (20,30] and 
(30,^). There are weak trends in the 3.1 and 3.8 /um 
bands to diminish in strength with greater optical depth. 
Symbols are as in Fig. 7. 
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Fig. 18. Change in slope around 3.2 and 3.5 /xm versus 
mass-loss rate. There is a trend for the spectral shape in 
the 3.1 3.7 /im region to turn from "bowl" to "cup", al- 
though the two MSX-selected carbon stars with the high- 
est mass-loss rates exhibit a "bowl" shaped spectrum. 
Symbols are as in Fig. 7. 
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Fig. 17. Equivalent width of the 3.1 jim band versus 
mass-loss rate before (left) and after (right) correction for 
dust emission, normalised to the peak in the distribution. 
Removing the dust emission creates a flatter and tighter 
distribution around unity for about half of the sample, 
but it creates large scatter for the remainder — many of 
which have low mass-loss rates. Symbols are as in Fig. 7. 



tion may in fact be seen against the dust continuum. If 
true this would imply a dense molecular atmosphere with 
a scale height of order a few stellar radii, which is seen in 
absorption against the dust emission arising from behind 
and adjacent to the limb of the star. The same may be 
true for the 3.8 fim band in some stars. 



A large molecular zone is also implied by the low tem- 
perature at the dust formation radius that was required 
to fit the SEDs. The molecular atmosphere must extend 
to at least the dust formation radius for dust to form, 
and a lower temperature at that radius implies that the 
dust forms farther from the star. If the temperature had 
been underestimated, then the veiling of the absorption 
bands by the dust continuum around 3-4 /xm would also 
have been underestimated. This would make the case even 
more convincing for part of the molecular absorption to 
be seen against background dust continuum emission. 

A weak trend can be seen in the differential spec- 
tral slope at 3.5 ^m as a function of mass- loss rate (Fig. 
18), suggesting that stars enduring higher mass-loss rates 
have a "cup" -shaped 3.2 -3.5 fim spectrum. At high mass- 
loss rate the 3.1 /xm band may on average be formed 
higher above the stellar photosphere, at a lower excita- 
tion temperature, which would explain the narrower shape 
("cup"). On the other hand, the two MSX-selected stars 
with the highest mass-loss rates defy such relationship. 

4.5. Pulsation and excitation conditions 

Pulsation is known to affect the absorption bands. The 
difficulty with investigating dependencies on pulsational 
properties is that the effects have a highly time-variable 
component. But the pulsation may also affect time- 
averaged properties, that we can study here. In addition 
to the stars with known pulsation properties listed in 
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• = 3.1 /xm 

o = 3.8 /xm. 
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Fig. 19. Equivalent widths at 3.1 (solid symbols) and 
3.8 fim (open symbols) versus Kg-band pulsation ampli- 
tude (left panel), and versus pulsational energy injection 
rate per cycle, E (right panel), in logarithmic units of Lq. 
Whilst no correlation is seen with amplitude, the equiva- 
lent width is larger for larger energy injection rates. 

Table 2, pulsation periods and amplitudes are known for 
IRAS 04496-6958 and IRAS 04557-6753 (Whitelock et 
al. 2003) and rough estimates for LI-LMC 1813 (van Loon 
et al. 2003). The L'-band amplitude of 1.6 mag for LI- 
LMC 1813 is scaled to a K^-band amplitude of AKg ~ 2.0 
mag (cf. Whitelock et al. 2003). 

The equivalent width (at 3.1 and 3.8 /xm) shows no 
correlation with the photometric amplitude, AKg (Fig. 
19, left panel; p = —0.15 and p = —0.21, respectively). 
However, the photometric amplitude is relative to the 
star's luminosity. An absolute quantity that more directly 
measures the amount of energy involved in the pulsation 
is the energy injection rate (van Loon 2002): 



L r exp (AKs/2.5) - 1 
2 ^ \ exp (AKs/2.5) + 1 



(1) 



This is the fraction of the bolometric luminosity that is 
used in expanding the star, assuming that the amplitude 
at Kg measures the bolometric amplitude and that the 
luminosity varies sinusoidally. It is by definition limited 
to < < L/2. The equivalent widths at 3.1 and 3.8 //,m 
show a weak correlation with E (Fig. 19, right panel; p = 
0.50 and p = 0.41, respectively), suggesting a positive link 
between stronger pulsation and a larger column density of 
the molecular atmosphere — in particular of C2II2. 

Pulsation also affects the excitation conditions. The 
fact that the 3.1 /zm band is narrower at high mass- loss 
rates is because both are related to the increased scale 
height as a consequence of strong pulsation. Indeed (Fig. 
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Fig. 20. Change in slope at 3.2 /xm (solid symbols) and 3.5 
/im (open symbols) versus Kg-band pulsation amplitude 
(left panel), and versus pulsation period (right panel). 
The 3.2-3.5 /im absorption wing is weaker for larger am- 
plitudes and longer periods, indicating a lower excitation 
temperature in a more extended molecular atmosphere. 

20), the differential spectral slope at 3.2 fim is weakly 
correlated with the Kg-band amplitude {p = 0.26) and 
pulsation period {p — 0.53) whilst the differential spectral 
slope at 3.5 /xm shows the opposite behaviour {p = —0.61 
and p = —0.45, respectively). This suggests that the 3.1 
lim absorption profile is narrower at larger photometric 
amplitude and longer pulsation period, supporting evi- 
dence for a lower excitation temperature in a more ex- 
tended molecular atmosphere. There is a possible com- 
peting mechanism associated with the pulsation: shocks 
induced by strong pulsation cause the band to broaden 
again (Yamamura et al. 1997). This might explain the 
"bowl" -shaped absorption in the two MSX-selected stars 
with the highest mass-loss rates. 

5. Discussion 

5.1. The population of dust-enshrouded carbon stars 

The population of visually bright carbon stars in the LMC, 
selected on the basis of objective prism surveys, is esti- 
mated to be A^vis - 9,000 (Kontizas et al. 2001). In the 
galaxy, many carbon stars become extremely obscured by 
the very opaque carbon dust they produce when they enter 
the superwind stage with M > 10~^ Mq yr~^ (e.g., Jura 
& Kleinmann 1990). These objects would escape detection 
in the optical and (shallow, wide-field) near-IR surveys 
in the LMC but could stand out at mid-IR wavelengths 
where their dust envelopes shine brightly. The IRAS Point 
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Source Catalogue was searched for obscured AGB stars in 
the LMC, and although several hundred such candidates 
were found (Loup ct al. 1997) only a few dozen were suc- 
cessfully identified with carbon stars, with another few 
dozen oxygen-rich massive AGB stars or red supergiants 
(Reid et al. 1990; Wood ct al. 1992; Zijlstra ct al. 1996; van 
Loon et al. 1997). This was but a small fraction of the ex- 
pected population of superwind AGB stars. However, the 
relatively deep near-IR imaging of IRAS positions by van 
Loon et al. (1997) recovered counterparts all the way down 
to their detection limit of Kg < 18 mag, whilst leaving half 
of their IRAS sample unidentified. Most of the faint coun- 
terparts had colours of carbon stars — later spectroscopic 
data confirmed their carbon star nature (van Loon et al. 
1999a; Trams ct al. 1999). Van Loon et al. (1997) took 
this as evidence for a hitherto uncovered population of 
very obscured carbon stars. 

The MSX Point Source Catalogue at 8.3 /v,m was cross- 
identified with the 2MASS JHKg catalogue to reveal the 
largest selection to date of candidate superwind AGB stars 
in the LMC (Egan ct al. 2001; Wood & Cohen 2001). 
Several hundred were found, with more than a hundred re- 
maining without near-IR counterpart. Given the result of 
van Loon ct al. (1997) and the shallower depth of 2MASS 
at Kg < 15 mag we can expect many of the MSX-selected 
candidate AGB stars and unidentified stellar MSX sources 
to be dust-enshrouded (IR) carbon stars. We here provide 
supporting evidence for this hypothesis, as our selection 
of targets from the MSX-2MASS compilation without op- 
tical counterparts turns out to almost entirely consist of 
IR carbon stars (at least 28 out of 30). 

Allowing for incompleteness in the MSX catalogue due 
to the detection threshold and confusion in crowded areas 
or areas with a complex background, the total population 
of IR carbon stars in the LMC may approach ~ 1, 000. 
This would suggest that, typically, a 1.5 3 Mq star spends 
around 10 per cent of its thermal-pulsing AGB lifetime be- 
ing heavily dust-enshrouded. With mass-loss rates in this 
superwind stage enhanced by a factor of 10 (or more) over 
those experienced during the unobscured stages, the mass 
return during the superwind stage of a carbon star must 
be (at least) 50 per cent of the total mass returned over 
its lifetime. This is in excellent agreement with the results 
found for cluster IR objects (van Loon et al. 2005b). 

5.2. Spectral diagnostics 

The 3.1 /um band is due to a mixture of HCN and C2H2, 
but it is difficult to determine the relative contributions 
of these molecules. In all MSX-selected IR carbon stars 
in the LMC the 3.8 yum band peaks at 3.75 /xm and is 
therefore probably due to C2H2; it never peaks at 3.9 ^m 
as it often does in the galaxy due to absorption by HCN 
and/or CS. The lower metallicity in the LMC would yield 
lower nitrogen and sulphur abundances, reducing the po- 
tential for forming HCN or CS. The equivalent widths of 
the 3.1 and 3.8 ^.m bands correlate very well, although the 



3.8 /xm band is not always seen whilst the 3.1 /xm is al- 
ways prominent. This suggests that when both bands are 
strong, the C2H2 molecule is abundant, whilst HCN may 
be more abundant when the 3.8 /xm band is absent. This 
is corroborated by the observation that the 3.57 /xm band 
due to HCN is seen both in objects with very strong 3.8 
/tm absorption and in objects without that C2H2 band. 
The 3.57 /tm shows no dependence on luminosity or mass- 
loss rate, suggesting that its abundance is primarily set by 
the initially available nitrogen with little evidence for ad- 
ditional nitrogen produced by the CNO cycle in more mas- 
sive AGB stars (those would probably not become carbon 
stars anyway due to Hot Bottom Burning). The strength 
of the 3.8 /tm band can thus be taken as a measure for 
the carbon overabundance resulting from the carbon pro- 
duced in the star itself, whilst the 3.57 /tm band measures 
the (initial) nitrogen abundance. 

The 3.1 /xm band sometimes shows a strong red wing 
of absorption extending up to 3.5 /tm. Although this could 
be due to the HNC isomer (Harris et al. 2003), it is not 
expected to form in great quantities under chemical equi- 
librium conditions. Unless stellar pulsation causes depar- 
tures from chemical equilibrium that result in enhanced 
HNC abundances, it is more likely that the broad absorp- 
tion is due to a high excitation temperature, closer to the 
stellar surface. The excitation conditions almost certainly 
depend on the strength of pulsation, and may vary as a 
function of phase in the pulsation cycle. We find evidence 
for the excitation temperature to be on average lower in 
stars with longer periods and/or larger amplitudes, indica- 
tive of a more extended molecular atmosphere. 

We have investigated in some detail the nature of the 
relatively sharp absorptions at 3.70-3.78 /tm. Of these, 
the 3.70 /xm "band" is the easiest to measure reliably. 
The 3.70 /tm and associated sharp features are seldomly 
as prominent in the SWS spectra of galactic carbon stars, 
suggesting that they may be indicative of the C enhance- 
ment seen in LMC carbon stars. The 3.70 /xm band is 
generally absent at low mass-loss rates whilst the 3.57 
/xm band is often quite strong in those stars, suggesting 
that the 3.70 /xm absorption is not due to HCN. Both 
CH and C2H2 have features in this spectral region. No 
CH absorption — either around 3.75 /xm or elsewhere in 
the spectrum — is evident in any of the IR carbon stars 
in the LMC. It is seen around 3.4 /tm in the optically 
bright carbon star MSX LMC 349 B (Fig. 4) and in rela- 
tively blue carbon stars in SgrD (Matsuura et al. 2005), 
but in neither case does absorption at 3.8 /tm occur. This 
suggests that in optically bright carbon stars CH domi- 
nates over C2H2 whilst in IR carbon stars the situation 
is reversed — probably a consequence of the higher C 
overabundance and/or denser atmosphere promoting the 
formation of larger carbon-chains. The 3.70 /xm band is 
strongest at intermediate-strong 3.8 /xm absorption, and 
these bands show an opposite behaviour with respect to 
the 3.2-3.5 /tm absorption. We suggest therefore that, like 
the broad 3.8 /xm band, the 3.70-3.78 /tm absorption com- 
ponents are due to C2H2, with excitation conditions de- 
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termining their relative strengths in the sense of the sharp 
features dominating at lower excitation temperature in an 
extended atmosphere and the broad absorption dominat- 
ing at higher excitation temperature in lower (or shocked) 
layers of the atmosphere. 

5.3. The assembly of larger molecules and dust grains 

The properties of the molcciilar atmosphere amongst the 
IR carbon stars in the LMC are very homogeneous, and in 
particular show little — if any — variation with luminos- 
ity or mass-loss rate. For instance, although C2H2 seems 
to be abundant amongst LMC carbon stars experiencing 
the superwind, once in that final evolutionary stage we 
were unable to identify a strong correlation between the 
C2H2 absorption and the mass-loss rate. The subsequent 
production of PAH molecules and carbon grains may be 
largely constant. Veiling of the 3-4 fim spectrum by the 
continuum emission of warm dust is shown to be impor- 
tant, but it is difficult to accurately account for. 

Strong pulsation and intense mass loss do go hand in 
hand, and we presented some evidence for an extended 
molecular atmosphere: low excitation temperature, satu- 
rated molecular bands, and a possible contribution from 
absorption seen against the dust envelope. Although this 
could be thought of as a physical link between the action of 
pulsation and the dust-driven wind, no clear cause-effect 
relationship is found. This might indicate that the pul- 
sation and atmospheric structure have already reached 
the threshold required for complete condensation into 
grains, leading to a "saturation" of the mass-loss mech- 
anism. That could explain the metallicity-insensitivity of 
the maximum mass-loss rate attained by AGB stars (van 
Loon 2000, 2002). 

Further insight into the molecule and dust formation 
should come from a more extensive comparison between 
oxygen-rich AGB stars and carbon stars in the LMC with 
those in the SMC, and between the properties of the 
molecular atmosphere seen at optical and infrared wave- 
lengths with the properties of the dust grains seen at mid- 
IR wavelengths (Cohen 1984). 

6. Summary of conclusions 

We presented 2.9-4.1 /zm spectra, obtained with ISAAC at 
the ESO/VLT, of a sample of 30 optically invisible stars 
selected from the MSX-2MASS compilation of Wood & 
Cohen (2001). At least 28 (over 90 per cent) turn out 
to be carbon stars. This confirms earlier indications (van 
Loon et al. 1997) for an extensive population of IR carbon 
stars experiencing the superwind stage at the end of the 
thermal-pulsing AGB. This population is now estimated 
to account for about 10 per cent of the thermal-pulsing 
AGB lifetime of LMC carbon stars, in agreement with 
findings from cluster IR objects (van Loon et al. 2005b). 

We presented an analysis of the molecular band 
strengths in the spectra of these and several other LMC 



carbon stars, and conclude that in the most evolved ob- 
jects the C2H2 abundance is large compared to that of 
HCN and CH. We interpret this as evidence for a low 
nitrogen abundance and a high carbon abundance as a re- 
sult of the low initial metallicity and efHcient production 
of carbon in LMC carbon stars, confirming and extending 
earlier such evidence (van Loon et al. 1999; Matsuura et 
al. 2002, 2005). 

No strong relationships are foimd between the strength 
of the molecular bands and either luminosity or mass-loss 
rate, but we found hints for the molecular atmosphere 
to be more extended at the highest mass-loss rates and 
strongest pulsation, characterised by a low excitation tem- 
perature and high C2H2 abundance. The pulsation, molec- 
ular atmosphere and dust wind of these IR carbon stars 
may have reached their maximum potential for driving 
mass loss. 
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